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Abstract 

We compute nonequilibrium dynamics of plasma instabilities in classical- 
Q\ i Statistical lattice gauge theory in 3+1 dimensions. The simulations are done 

O I for the first time for the SU(3) gauge group relevant for quantum chromo- 

^ ■ dynamics. We find a qualitatively similar behavior as compared to earlier 

investigations in SU(2) gauge theory. The characteristic growth rates are 
\ about 25% lower for given energy density, such that the isotropization pro- 

cess is slower. Measured in units of the characteristic screening mass, the 
primary growth rate is independent of the number of colors. 

1 Introduction 

(D ■ The theoretical understanding of the apparent fast thermalization in collision ex- 

r^, ' periments of heavy nuclei at the Relativistic Heavy Ion Collider provides a chal- 
lenge for theory. It was noticed that plasma instabilities in the locally anisotropic 

^ ! medium might play an important role for a rapid isotropization of the equation of 

1^ I state, which is relevant to explain the observed hydrodynamic behavior [[II 121 El 

^ I HI [3. Extensive studies have been carried out for SU(2) pure gauge theory using 

Cf) • the hard-loop effective theory of soft excitations, which is based^on coUisionless 



kinetic theory for hard particles coupled to a soft classical fieldll] This approach 
^ . neglects quantum corrections and may also be considered as an approximation 

I of the classical-statistical field theory limit of the respective quantum gauge the- 

ory [[UHl. Classical-statistical lattice gauge theory provides a quantitative descrip- 
tion in the presence of sufficiently large energy density or occupation numbers per 
mode. The simulations are done by numerical integration of the classical lattice 
equations of motion and Monte Carlo sampling of initial conditions. 

In this work we present for the first time classical-statistical lattice gauge the- 
ory simulations for the SU(3) gauge group in 3-1-1 dimensions, relevant for quan- 
tum chromodynamics (QCD). Simulations using the hard-loop approximation for 
SU(3) gauge group in l-i-l dimensions have been done in Ref. [[51 . Though the 
dynamics in one and three spatial dimensions are known to be very different, it 



X 



'For numerical simulations which take into account the backreaction of the soft fields on the 
hard particles using a Boltzmann-Vlasov treatment see Studies using transport or kinetic 
equations were also carried out |[6l in the spirit of the earlier bottom-up scenario Q, where no 
instability is present. 
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is not expected that SU(3) results will qualitatively change as compared to those 
previously obtained for SU(2), however a quantitative estimate of the involved 
time scales seems imperative. This work is a follow-up to our previous studies 
concerning the numerically less demanding gauge group SU(2) ||4l. In this Let- 
ter we only describe the relevant changes and refer to that reference for further 
computational details. We find a qualitatively similar behavior as compared to the 
earlier investigations in SU(2) gauge theory, but for given initial energy density 
the characteristic growth rates are about 25% lower such that the isotropization 
process is slower. 

The paper is organized as follows. In Section 2, we outline the setup of our 
calculations, and the algorithm for solving the equations of motions for the SU(3) 
group. In Section 3 we present the results of the calculations and explain the 
dependence of growth rates on the gauge group in terms of diagrammatic s. We 
conclude in Section 4. 



2 Classical-statistical gauge field theory on a lattice 

Following Ref. dH we use the Wilsonian lattice action on a 3-1-1 dimensional 
Minkowskian lattice, 

S[U] = -^oJ]J]{^(Trf/,o,+Trf/g-l| 

+^^ZZ{l^(T^^-V + Trf/:,.)-l|, (1) 

written in terms of plaquette variables 

Ux,,y^U,,,U,,l,,MU^Ul, (2) 

where f/J = Ux^^^y . Here [7^,^/ denotes the link variable which is the parallel 
transporter associated with the link from the neighboring lattice point x + fi to 
the point x = (t,x) in the direction of the lattice axis ft = 0, 1,2, 3. The lattice 
parameters are defined as 

27Trl 2Trl 

^0 Sir 

where y = ajat is the ratio of the spatial and temporal lattice spacings, and we 
will consider = gs = g as the coupling constant of the lattice theory. 
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Varying the action ([T]) w. r. t. the spatial link variables U^^j yields the leapfrog- 
type equations of motion 

(4) 

where /l^ are the generators of SU(A'^) satisfying TriA^A^) = ^5"^. In the case of 
SU(2) these can be related to the Pauli matrices cr" through A" = ^cr". For SU(3), 
the generators are given by the Gell-Mann matrices divided by two. The electric 
fields in Eq. (H)) are 

E'(x) = — ImTr (t/,,0,/) . (5) 
Varying the action ^ w. r. t. to a temporal link gives the Gauss constraint 

j][E^(x)-Ul..E';(x-j)U,_jj]=0. (6) 

Using the gauge freedom the dynamics is computed in temporal axial gauge, 
which is defined by the gauge condition i7v,o = 1- This choice still leaves the 
possibility of a time-independent gauge fixing which we do not specify. The al- 
gorithm for the numerical simulation may be summarized as follows: @ 

1. Using E^.{t-ai, x) and t/(/,x),j the electric field is evolved to x) with dH). 

2. The temporal plaquette t/(?,x),oy 6 SU(A'^) is calculated which satisfies 
ImTr([/(,,,),o;^') = {asatgE){t,x). 

3. The link variable U(t+a,.x).j is determined at the next time step from the defin- 
ing equation of the temporal plaquette in temporal axial gauge: U(t+a,,x),j = 

U(t,x),OjU{t,x),j- 

The main difference between the algorithms for SU(2) and SU(3) occurs in step 
2. For U 6 SU(2) and given ImTr(i7cr") = 2ba, with a = 1, 2, 3, one can represent 
U as U = VI - babat + ibacr" for the solution which is close to the unit matrix. 

For SU(3) (and SU(A^) with > 3 in general) there is no simple corresponding 
procedure in step 2 and we solve for the temporal plaquette numerically. For a 

-Another simulation algorithm would be the Hamiltonian algorithm pi, which keeps track of 
the electric fields E'j(x) in the Lie algebra space instead of the link variables in group space. This 
way one would need to calculate exponentials of matrices to calculate the time evolution of the 
spatial links. 
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general complex 3x3 matrix, which we use in the program, there are 18 real 
variables (3x3 complex numbers) and 28 real, but not independent, equations 
(8 from ([5]), 2 from detUxflj = 1, and 18 from UxfljUl^j = 1). Since one only 
needs 18 equations, we use part of the unitarity equation: components 12, 13, and 
23, and the real part of components 1 1 and 22 because the imaginary parts do not 
give any constraints on the variables. These equations are solved with the multi- 
dimensional Newton method |[8l . The fastest simulation speed was achieved by 
using the unit matrix as a starting point of the Newton method, which typically 
leads to convergence in about 3-4 steps. Since for every new link variable one has 
to iterate the Newton method in 18 dimensions a few times, the SU(3) simulation 
is considerably slower than the corresponding SU(2) simulation. 
The initial Gaussian probability functional is chosen as in Ref. IIH such that 

~2 2 2 2 

where the gauge fields are calculated from the link variables using 

A'it, X) = — ImTr(f/(,,),,a^). (8) 

We consider E" = -A" = at f = fulfilling the Gauss constraint Q. We typ- 
ically choose A- «; A, and the distribution is practically 6{p^)-like on the lattice. 
Here A determines the typical transverse momentum of the gluons and may be as- 
sociated with the saturation scale Qs at time 27 ' in the saturation scenario We 
solve the equations of motion with the aforementioned algorithm for a set of ini- 
tial configurations sampled according to Eq. (|7]) and compute expectation values 
as averages over the results from the individual runs. 
The local energy density in lattice units is determined by the action as 

2 Tr [/.,,oy + Tr U^.jo ^ Tr U_^jk + Tr U^^kj ^ 
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sit,x) = — 



^ 6 ^2^^ 



(9) 



and we denote its average value hy s = {s{t,x) ). We choose the factor A ap- 
pearing in Eq. ^ such that s = 0.05. For the conversion to physical units we 
fix the lattice spacing as from the relation between the physical average energy 
density and its lattice analogue according to e = £-a~'^. If g is taken to be different 
from one then the lattice spacing is altered by a factor Ij yfg , which follows from 
Eq. dll). The values for will not be altered significantly as long as g ~ 0{\). 
For later reference, we note that the relation between s, N and A is approximately 
given by 

4 N--\ 

s(t,x) « y Z {(^i^^''^))' + - (N' - (10) 
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Figure 1: Fourier coefficients of ttie squared modufus of ttie gauge fiefd versus 
time for tliree different momenta paraiiei to tlie z-axis. Compared is the time 
evoiution of SU(2) and SU(3) gauge fieids for tlie same energy density e. 



for our initiai conditions, wliere tlie magnetic fields are calculated from the spatial 
plaquettes according to 

B){x) = £^'''4- ImTr (t/,,,//) , dD 
similarly to the electric fields dH). 

3 Results 

Fig. [T] shows the nonequilibrium time evolution of the color- averaged squared 
modulus of two different Fourier coefficients A{t, p) of the gauge field in three 
spatial dimensions. They are displayed as a function of time, normalized by the 
corresponding field values at initial time, where time is measured in appropriate 
units of the initial energy density e. The lattice size is 64^ and the initial transverse 
width A = 1.06 e^^"^. Here { |A(?, p)p ) may be associated to a particle number 
divided by frequency. 

For comparison we show the SU(3) results together with the corresponding 
SU(2) results. The behavior of both gauge groups is qualitatively very similar. 
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Figure 2: The primary and secondary growth rates for the SU(2) and SU(3) gauge 
theory using the same energy density. The employed lattice sizes are 64^ and 96^. 
The exponential growth was fitted from an average of ca. 20 and ca. 100 runs for 
the SU(3) and SU(2) gauge groups, respectively. 



The plotted low-momentum modes clearly show exponential growth starting at 
the very beginning of the simulation. In contrast to these "primary" instabilities 
occurring at low momenta, one observes from Fig.[I]that gauge field modes at suf- 
ficiently high momenta do not grow initially. The higher wave number modes typ- 
ically exhibit exponential growth at a "secondary" stage that sets in later, but with 
a significantly larger growth rate. The secondaries arise from fluctuation effects 
induced by the growth in the lower momentum modes, which can be explained by 
taking into account (2PI) resummed loop diagrams beyond the hard-loop approxi- 
mation as discussed in detail for SU(2) gauge theory in Ref. IIH. In that case it has 
been shown that the exponential growth saturates when all loop diagrams become 
of order one, which leads to a subsequent slow evolution towards a power-law 
regimeEl Since the fastest growing mode in SU(2) gauge theory grows earlier to 
sizes where loop effects for higher modes become important, one observes from 
Fig.[I]that the secondary growths also start earlier compared to the SU(3) case. 
Fig. [2] displays the momentum dependence of the growth rates for ( \A(t, p)p ) 

^For related discussions in the context of scalar inflaton dynamics in the early universe see 
Refs. HOllIT] and Cni for gauge fields. 
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Figure 3: The tadpole diagram whose non-vacuum part determines mr as de- 
scribed in the main text. 



^SU(3) ~ Q ^SU(2) ■ (13) 



obtained from a fit to an exponential, which is done separately for the primary and 
secondary growth rates. One observes that while the primary rates are approxi- 
mately 25 % bigger for SU(2), the secondary rates can be similar within the given 
errors. The dependence of the primary growth rates on the gauge group can be 
understood from a diagrammatic analysis as follows. As outlined in Section [2l 
the overall gauge field amplitude A defined in ([7]) is chosen such that the aver- 
age lattice energy density s has a prescribed value. For simulations of difi'erent 
gauge groups SU(Ni) and SU(A^2) at the same s using initial conditions sampled 
from dV]) with identical A and A-, Eq. (fTUl) implies that the respective overall gauge 
field amplitudes A are related according to 

(A^i-1)AL(;vo = (A^2-1)AL(;v,)- (12) 

Plugging in the relevant numbers we find for SU(2) and SU(3) 

3 
8 

In Ref. [4J we have analyzed for SU(2) gauge theory relevant diagrammatic 
contributions to the nonequilibrium evolution at various characteristic time scales, 
which holds along the same lines for SU(3). As a characteristic self-energy con- 
tribution relevant for the primary growth rates at early times, we consider the 
non- vacuum part of the tadpole diagram shown in Fig. [31 For SIJ{N) gauge 
theory the parametric A^^-dependence of that contribution is 

'^T,SV(N) ~ '^^SU(A') ' (14) 

which implies 

^r,su(Wi) _ Ni ^su(M) 

2 - ju n ^^^> 

"*r,su(7V2) 2 ^su(iV2) 
at fixed e. Using Eq. (fT3l) we can thus compare the change in mj when going from 
SU(2) to SU(3) at fixed energy density which amounts to 

3 

/"T.SUO) = ^'«T,SU(2)- (16) 
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Figure 4: Comparison of classical-statistical simulation results for the primary 
growth rates as a function of momentum for SU(2) and SU(3) gauge theories, 
measured in units of the respective my. The same parameters are used as in Fig.[IJ 



Fig. m shows that the classical-statistical simulations for SU(2) and SU(3) gauge 
theory give indeed approximately the same results for primary growth rates in 
units of the respective 

In order to obtain an estimate in physical units one may consider an initial en- 
ergy density of about 5-25 GeV/fm^ for RHIC experiments, and a projected factor 
of about two more for LHC energies. The inverse of the maximum primary growth 
rate for | A{t, p) |^ then leads to the characteristic time scales 

r;;L.pr. - 1.6 - 2.4fm/c (RHIC), (17) 
r;;lx.pr. - 1-3 - 2.0fm/c (LHC). (18) 

The results are rather insensitive to the precise value of the initial energy den- 
sity because they scale with the fourth root of e. For comparison, a time scale 
associated with the largest observed secondary growth rate is about a factor of 

"^This mass may be directly related to nica in the 'hard loop' approximation, calculated from the 
anisotropic distribution of the hard modes ( 1 1 . Note that for our case nir does not scale with like 
the thermal mass, because we keep the energy density fixed, which results in different temperatures 
for theories with different numbers of degrees of freedom. 
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three shorter than what is given in (flTl) . However, even though secondaries can 
reach considerably higher growth rates than primaries, they start later. As a conse- 
quence, a certain range of higher momentum modes can 'catch up' with initially 
faster growing infrared modes before the exponential growth stops, as seen in 
Fig.[il This leads to a relatively fast effective isotropization of a finite momentum 
range, while higher momentum modes do not isotropize on a time scale charac- 
terized by plasma instabilities. 

We study the process of isotropization in terms of the local energy density ^ 
which is a gauge invariant quantity. For this, we Fourier transform s{t, x) with 
respect to x and compute the absolute value of the ratio 

with PL II z , Pt- ± z , and I Pl I = I pr I . (19) 

If at some time t the system reaches an isotropic state the mean absolute value 
of (fT9l ) has to be one for all momenta. The time evolution of (fT9l ) is depicted in 
Fig.[5]for several momenta. Similarly to the SU(2) case one observes that the low- 
momentum modes of the energy density isotropize, while the high-momentum 
modes are still anisotropic after the instability growth ends. 

Inspecting Fig. |5] one sees that the Fourier coefficients of s become isotropic up 
to approximately | p | - 1-2 g^^"* at the time of saturation. Using the same values 
for the physical energy density s as for (fTTl) one finds efl'ective isotropization up 
to a characteristic momentum of about 

Ipl < IGeV, (20) 

which approximately agrees with the case of the SU(2) gauge theory. 



4 Conclusions 

Nonabelian plasma instabilities in SU(3) gauge theory relevant for QCD exhibit 
a qualitatively similar behavior as previously observed for the SU(2) group. The 
main quantitative differences concern the reduction of primary growth rates by 
about 25 % for given initial energy density. These differences can be related to 
the parametric dependence on the number of colors of one-loop corrections to the 
self-energy. As a consequence, we find that the properly rescaled primary growth 
rates agree even quantitatively rather well for the different gauge groups. Though 
the nonlinear dynamics underlying the secondary growth rates leads to remain- 
ing quantitative differences between SU(3) and SU(2), this fact is less important 
for phenomenology. The primary growth rates determine the characteristic time 
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Figure 5: The ratio of the longitudinal and transverse modes of the Fourier trans- 
form of the spatially dependent energy density for various momenta, measured on 
= 64^ lattices, using SU(3) gauge group. 



scales for isotropization at low momenta, which turns out to be (even) slower than 
previously suggested by SU(2) results. 
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